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PROCEDE POUR EVALUER UNE GRANDEUR PHYSIQUE REPRESENTATIVE D'UNE INTERACTION ENTRE 
UNE ONDE ET UN OBSTACLE. 

L'invention concerne la moderation des interactions 
entre une onde incidente et un obstacle, notarnment dans le 
domaine du controle non destructif. Selon I'invention, on 
maille la surface de I'obstacle et Ton attribue au moins une 
source (Si) a chaque element de surface (dS f ). On calcule 
ensuite des conditions aux limites en chaque maille de 
('obstacle et Ton en deduit des valeurs de source. A partir 
d'une matrice d' interaction et de ces valeurs de sources, on 
estime une grandeur physique representative de I 1 interac- 
tion entre Ponde et I'obstacle, en tout point de I'espace. 
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Procede pour evaluer une grandeur physique rep resentative 
dMane interaction entre une onde et un obstacle 

L 1 invention concerne ia mod61isation des interactions 
5 entre une onde incidente et un obstacle de cette onde, 
notarament dans le domaine du contrdle non destructif. 

La mod61isation des interactions entre une.. onde et un 
obstacle recevant cette onde, tel qu'une cible plac£e dans 
10 la zone sensible d'un capteur, trouve une application 
avantageuse en contrdle non destructif. 

On connait un procede de modfelisation dit "par 616ments 
finis" consistant a appliquer un pavage de l f espace 
15 tridimensionnel entourant l 1 obstacle et a evaluer les 
interactions precipes pour tous les pav6s de l'espace. 

Les procedes de calcul par "Elements finis" apportent une 
resolution a un probldme pose sous forme d' Equations 

20 diff§rentielles partielles. lis se fondent sur une 
representation de l'espace d' 6tude par un assemblage 
d' Elements finis, a l'interieur desquels des fonctions 
d' approximations sont d§terminees en termes de valeurs 
nodales de la grandeur physique recherch6e. Le probleme 

25 physique continu devient done un probleme discret aux 
elements finis oix les grandeurs nodales sont les nouvelles 
inconnues- De tels proc6des cherchent done a approcher la 
solution globale, plutdt que les equations aux d£riv6es 
spatiales partielles de depart. 

30 
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La discretisation de l'espace pris en compte assure que ce 
dernier soit entitlement recouvert par des elements finis 
(lignes, surfaces ou volumes), cette operation porte le 
nom de "maillage" dans l'espace bidimensionnel (2D) ou de 
5 "pavage" dans l'espace tridimensionnel (3D). Les elements 
mis en jeu sont soit rectangulaires ou triangulaires en 
2D, soit parall£16pip6diques ou tetraedriques en 3D. lis 
peuvent etre de tailles difftrentes, repartis unif ornament 
sur la surface ou non. 

i° 

En g£n6ral, la grandeur physique recherche, telle qu'un 
potentiel electrostatique ou une valeur de pression, est 
connue sur la' front id re du domaine. Cette frontiere peut 
etre fictive. On y impose des conditions aux limites. Le 
15 potentiel est done inconnu a 1' interieur du meme domaine. 
On d£finit alors un noeud comme etant un sommet d' un 
Element. Les inconnues du probleme sont done les valeurs 
du potentiel en chaque noeud de 1' ensemble du domaine. 

20 A titre d' illustration, la figure 6 de I'art anterieur 
reprtsente un exemple de surface, constitute par deux 
materiaux Ml et M2, de proprietes 61ectromagnetiques 
diff£rentes, et maillte par des elements triangulaires 
comportant chacun trois nceuds Ai, Bi et Ci. L 1 ensemble du 

25 domaine est delimits par une frontiere F. 

Une fois le maillage dtfini, plusieurs approches existent 
pour transformer la formulation physique du probleme en 
une modelisation discrete aux elements finis. Si le 
30 probl&me est formula par des equations dif f trentielles et 
consiste & minimiser une f onctionnelle, on applique 
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generalement une methode variationnelle. Cette 
transformation aboutit £ une formulation matricielle dont 
la resolution donne les solutions nodales, les solutions 
aux points non nodaux etant obtenues par interpolation 
5 linfeaire. 

Neanmoins, de tels calculs, en trois dimensions , 
requierent des ressources inf ormatiques iraportantes et 
generent des temps de calcul tr£s longs , en d6pit de 
10 1 1 accroissement des performances de logiciels permettant 
d 1 implementer ces calculs. 

Certes, les probiemes 2D, souvent simplifies par des 
conditions de symetrie avantageuses pour ne mod61iser 

15 qu'une partie de la geom6trie, se resolvent rapidement. 
Mais il n'en est rien pour les probiemes 3D, qui sont les 
plus frequents. La figure 6 montre combien la finesse du 
maillage, c'est a dire le rapport entre la taille d' un 
element et celle du plus petit detail du domaine, fait 

20 croitre le nombre de noeuds. 

En consequence, la quantity d' equations et d'inconnues 
augmente proportionnellement , et, de la, le temps de 
calcul necessaire a la resolution du problfeme. II est 
25 important de preciser que la generation du maillage, a 
savoir la discretisation de l'espace de travail, et la 
generation de la liste des noeuds consomme un temps de 
calcul superieur a celui necessaire ci la resolution du 
probleme . 

30 

La presente invention vient ameiiorer la situation • 
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Elle propose a cet effet un proc6d§ pour evaluer une 
grandeur physique associee a une interaction entre une 
onde et un obstacle, dans une region de l'espace 
5 tridimensionnel, dans lequel : 

a) on determine, par maillage, une pluralite 
d f 6chantillons de surface dont une partie au moins 
reprfesente la surface d'un obstacle recevant une onde 
principale et £mettant, en r^ponse, une onde secondaire, 

10 et l'on attribue a chaque 6chantillon de surface au moins 
une source 6mettant une onde elementaire repr£sentant une 
contribution k ladite onde secondaire, 

b) on forme un systeme matriciel comportant : 

- une matrice d' interaction, inversible, appliqu&e a 
15 une region donn6e de l'espace et comportant un nombre 

de colonnes correspondant a un nombre total de 
sources, 

une premiere matrice colonne dont chaque coefficient 
est associ6 & une source et caract6rise l'onde 
20 c ' elementaire qu'elle 6met, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d f une grandeur physique representative de 

25 l'onde emise par 1' ensemble des sources en ladite 

region donn6e, 

c) pour estimer les coefficients de la premiere matrice 
colonne, on affecte des valeurs de grandeur physique 
choisies a des points predetermines, associ§s chacun k un 

30 echantillon de surface, lesdites valeurs choisies etant 
rang^es dans la seconde matrice colonne, et l'on multiplie 
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cette seconcle matrice colonne par l 1 inverse de la matrice 
d' interaction appliquee auxdits points predetermines, 
d) pour 6valuer ladite grandeur physique repr6sentant 
l'onde £mise par 1' ensemble des sources dans une region 
5 donnee de l'espace tridimensionnel, on applique la matrice 
d' interaction a ladite region donnee et on multiplie cette 
matrice d' interaction par la premiere matrice colonne 
comportant les coefficients estimes a l'etape c) . 

10 Ainsi, selon l'un des avantages que procure la presente 
invention, l'etape de maillage a) ne concerne qu'une ou 
plusieurs surfaces, tandis que le procede de modelisation 
de type "par elements finis" n6cessite un pavage de tout 
l'espace avoisinant 1' obstacle, ce qui permet de r6duire, 

15 dans la mise en oeuvre du procede selon 1* invention, les 
ressources memoires et les temps de calcul n^cessaires. 

Le proc6de selon 1 ' invention s ' applique aussi bien S une 
onde principale emise par une source lointaine qu'a une 
20 onde principale emise en champ proche. 

Avantageusement, pour ^valuer une grandeur physique 
representative d'une interaction entre un element 
rayonnant une onde principale en champ proche et un 

25 obstacle recevant cette onde principale, 

- a l'etape a), on determine en outre, par maillage, une 
pluralite d ' echantillons de surface representant ensemble 
une surface active de 1' element rayonnant l'onde 
principale, et l'on attribue a chaque echantillon de la 

30 surface active au moins une source emettant une onde 
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felementaire repr^sentant une contribution a ladite onde 
principale, 

- on applique les etapes b) , c) et d) aux echantillons de 
la surface active, et 

5 on lvalue ladite grandeur physique reprfesentant 

l f interaction entre l 1 element rayonnant et 1" obstacle dans 
une region donnee de l T espace tridimensionnel f en prenant 
en compte la contribution, en ladite region donn6e, de 
l'onde principale emise par 1' ensemble des sources de la 

10 surface active et la contribution de l'onde secondaire 
emise par 1' ensemble des sources de la surface de 
1 ' obstacle . 

On entend par les termes "element rayonnant" aussi bien un 
15 emetteur de l'onde principale, tel qu'un g§n£rateur 
d'ondes, qu'un rfecepteur de l'onde principale, tel qu'un 
capteur de cette onde. 

Dans, un premier mode de realisation, la grandeur physique 
20 a evaluer est une grandeur scalaire et, a l'§tape a), on 
attribue une source unique a chaque <§chantillon de 
surface. 

Dans un second mode de realisation, la grandeur physique a 
25 ^valuer est une grandeur vectorielle exprim§e par ses 
trois coordonn£es dans I'espace tridimensionnel, et on 
attribue, A 1'etape a), trois sources a chaque 6chantillon 
de surface. 

30 Dans une realisation avantageuse, pour estimer, a 1'etape 
d) , la contribution de l'onde secondaire dans la region 
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donn^e de l'espace, les valeurs de grandeur physique 
choisies a l'6tape c) sont fonction d'un coefficient 
predetermine de reflexion et/ou de transmission de I'onde 
principale par chaque echantillon de surface de 
5 1' obstacle. 

Ainsi, on comprendra que I'onde secondaire peut aussi bien 
correspondre £ une reflexion de I'onde principale, qu'a 
une transmission de I'onde principale , ou encore a une 
10 diffraction de I'onde principale. Dans cette realisation 
avantageuse, l'etape c) correspond finalement a une 
determination des conditions aux limites a la surface de 
1' obstacle, en tant qu' interface entre deux milieux 
distincts notamment dans une het^rostructure . 

15 

Par ailleurs, pour un contrdle non destructif d'une cible 
formant un obstacle d'une onde principale, on attribue un 
coefficient de reflexion ou de transmission choisi & tous 
les points predetermines de la surface de la cible, et on 

20 compare une simulation obtenue par la mise en ceuvre du 
procede au sens de 1' invention avec une mesure 
experimentale. Ainsi, les points de la surface de la cible 
qui, dans la mesure experimentale, ne verifient pas la 
simulation correspondent a des inhomogeneites ou a des 

25 impuretes sur la surface de la cible. 

Dans une autre approche, on connait les proprietes 
globales de 1' obstacle, notamment en transmission et/ou en 
reflexion. Par la mise en ceuvre du proced6 de 1' invention, 
30 on optimise alors la position dans l'espace d'un capteur 
ou m§me la forme de ce capteur, pour une application ci un 
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controle non destructif , ce capteur etant destine a 
analyser une cible forraant un obstacle de l'onde 
principale. 

5 A cet effet, dans une realisation avantageuse, on compare 
une plurality de valeurs de la grandeur physique estim6e a 
l'£tape d) du proc6de au sens de 1' invention, obtenues 
pour une pluralite de regions de l'espace, pour 
select ionner une region candidate pour la disposition d'un 

10 element rayonnant destine a interagir avec l 1 obstacle. 

Coitime indique ci-avant, on entend par les termes "61£ment 
rayonnant" aussi bien un capteur qu'un gen6rateur de 
l'onde. On comprendra ainsi que 1 ' optimisation de la 

15 position de 1' element rayonnant peut s'appliquer aussi a 
1' optimisation de la disposition ou de la forme d'un 
gen§rateur d'ondes. Par exemple, la presente invention 
trouve encore une application avantageuse a la disposition 
de haut-parleurs dans uri volume ferm£, d61imite par des 

20 obstacles, comme par exemple l'habitacle d'un vehicule 
automobile . 

D'autres caracteristiques et avantages de 1' invention 
apparaitront a l'examen de la description d£taill6e ci- 
25 apres, et des dessins annexes sur lesquels : 

- la figure 1A repr£sente sch6matiquernent les surfaces 
respect ives d'un element rayonnant ER 6mettant une onde 
et d'un obstacle OBS recevant cette onde, maillees en 
vue d'evaluer une grandeur scalaire representative de 
30 l'onde en un point M de 1'espace tridimensionnel ; 



la figure IB reprEsente en detail un Echantillon de 
surface dSi correspondant a une maille de la figure 1A, 
ainsi qu'une source Si associEe a 1' Echantillon de 
surface dSi ; 

la figure 2A reprEsente schEmatiquement les surfaces 
respect ives d'un ElEment rayonnant ER emettant une onde 
et d'un obstacle OBS recevant cette onde/ maillees en 
vue d'Evaluer une grandeur vectorielle representative 
de l'onde en un point M de l'espace tridimensionnel ; 
la figure 2B reprEsente en detail un Echantillon de 
surface dSi correspondant £ une maille de la figure 2A, 
ainsi que trois sources associees SAi, SBi et SCi ; 
la figure 2C represents, en vue de face, une surface 
maiilee dont chaque Echantillon de surface comporte 
trois sources SA±, SBi et SCi, pour 1 1 estimation d'une 
grandeur vectorielle ; 

la figure 3A represente, a titre illustratif, les 
armatures d'un condensateur, de potentiels Electriques 
respect if s VI et V2, pour 1' estimation d'un potentiel 
Electrique au point M de l'espace tridimensionnel, & 
chaque Echantillon de surface dSi de la figure 3A Etant 
associEe une unique source Si ; 

la figure 3B reprEsente, & titre illustratif, les 
armatures d'un condensateur, de potentiels Electriques 
respectifs VI et V2, pour l'estimation d'un champ 

Electrique E(M) , au point M de l'espace 
tridimensionnel, & chaque Echantillon de surface dS ± de 
la figure 3B Etant associEe trois sources SAi, SBi et 
SCi ; 
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- la figure 4A representee comme les figures 1A et 2A, 
une interaction entre un Element rayonnant ER et un 
obstacle OBS, pour evaluer une grandeur physique 
(scalaire ou vectorielle) en un point M dans une 
5 portion de l'espace delimitee par la surface de 

1' element rayonnant et la surface de 1' obstacle, ce 
point de l'espace M recevant a la fois l'onde 6mise par 
1' element rayonnant et l'onde r6flechie par 
1' obstacle ; 

10 - la figure 4B, complementaire de la figure 4A, 
repr^sente une transmission par 1' obstacle OBS de 
l'onde emise par l f element rayonnant ER, en un point M 
d'un demi-espace d£limit6 par le plan que forme la 
surface de 1' obstacle OBS ; 

15 - la figure 5A repr6sente schematiquement un obstacle 
OBS, de dimensions finies, avec des sources associ6es 
aux 6chantillons de surface agencies pour estimer une 
grandeur representative d'une reflexion de l'onde sur 
1 ' obstacle ; 

20 - la figure 5B, en complement de la figure 5 A, reprfesente 
schematiquement un obstacle OBS, de dimensions finies, 
avec les sources associees aux 6chantillons de surface 
dispos6es pour estimer une grandeur representative de 
la transmission de l'onde par 1' obstacle ; 

25 - la figure 5C repr6sente une simulation d'une onde 
ultrasonore 6mise par un Element rayonnant ER et se 
propageant vers un obstacle OBS ; et 
- la figure 6 repr6sente un maillage de milieux 
tridimensionnels, pour 1 ' application d'un procede de 

30 calcul par "Elements finis", au sens de l'£tat de la 

technique. 
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On se refere tout d'abord a la figure 1A, sur laquelle la 
surface d'un obstacle OBS, recevant une onde, est mailiee 
selon une pluralite d' echantillons de surface dSi a dS 4 , 
5 conform6ment a l'etape a) precit^e. 

En se referant a la figure IB, a chaque echantillon de 
surface dSi est associe un hemisphere HEMi, tangent ci 
1 1 echantillon de surface dSi en un point de contact Pi- 

10 Pref erentiellement, ce point de contact Pi correspond au 
sommet de 1 ' hemisphere HEMj.. Pour 1' estimation d T une 
grandeur physique scalaire au point M (tels qu'un 
potentiel electrostatique, une pression acoustique ou 
autre) , une source unique Si est associee a 1 ' echantillon 

15 de surface dSi. Corame on le verra plus loin, dans le cas 
de 1' estimation d'une grandeur vectorielle dans un point 
de l'espace M, on affectera plutot trois sources £ chaque 
echantillon de surface dSi, 

20 Preferentiellement, 1 hemisphere HEMi est construit comme 
decrit ci-apres. Pendant l'etape de maillage a) precitee, 
on evalue, d'une part, la surface de 1' obstacle OBS et, 
d' autre part, on choisit un nombre d T echantillons de 
surface dSi selon la precision souhaitee de 1' estimation 

25 de la grandeur physique au point M. Ainsi, la surface d'un 
echantillon dSi est donnee par S 0 /N, ou S c correspond £ la 
surface totale de l 1 obstacle et N correspond au nombre 
choisi d' echantillons de surface dSi. 
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L' hemisphere HEMi est de meme surface que 1 1 echantillon 
dSi. Ainsi, le rayon Ri de l'hemisph^re se deduit de 
l 1 expression : 

Chaque maille que represente un echantillon de surface dSi 
presente, dans l'exemple decrit, une forme de 
parallelogramme, de centre Pi correspondant au point 
d f intersection des diagonales de ce parallelogramme. 
L 1 hemisphere HEMi est tangent £ 1 1 Echantillon de surface 
dSi en ce point Pi. Bien entendu, les mailles peuvent etre 
de forme differente, triangulaire ou autre- On indique de 
fa<?on g6n6rale que le point Pi correspond au barycentre de 
la maille. 

On definit ainsi la position de la source Si (situ6e au 
centre de 1' hemisphere HEMi). La distance s^parant la 
source Si du point de contact Pi correspond au rayon R t de 
1 T hemisphere HEMi et la droite qui passe par les points Pi 
et Si est orthogonale & la maille dSi. 

Dans l'exemple represents sur la figure 1A, on maille en 
outre la surface d'un element rayonnant ER, correspondant 
par exemple a un g6n£rateur d'ondes. A chaque maille de la 
surface de l f element rayonnant est associe un 6chantillon 
de surface dS f if comme on le verra plus loin. 

Le systeme matriciel que l'on met en forme & I'etape b) 
precitee correspond a : 
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V(M 2 ) 


= F x 




,V(M N ), 







[1] 



ou : 

- les coefficients Vj (avec j=l,2,...,N) de la premiere 
matrice colonne correspondent a des valeurs associees 
chacune & une source, telle qu'une charge §lectrique 
(dans le cas de 1 1 estimation d ? un potentiel 
61ectrique) , ou encore a un flux magnetique (dans le 
cas de l'estimation d'un potentiel magnetique) , ou 
encore a une Vitesse sonore (dans le cas de 
1' estimation d'une pression acoustique li6e a la 
propagation d'une onde sonore) ; 

- les coefficients V(M±) (avec i=l,2,...,N) de la seconde 
matrice colonne correspondent chacun a une valeur de la 
grandeur physique (un potentiel 61ectrique ou 
magnetique ou une pression) £ estimer en un point M± de 
l'espace ; 

- la matrice d' interaction F comporte des coefficients 
C±,j dont 1' expression g§n6rale est donn6e par : 

Gj=f{M,Sj) [2] 

On comprendra ainsi que les coefficients de la matrice F 
sont des coefficients d 1 interaction qui dependent de la 
distance s6parant chaque point de l'espace Mi d'une source 
Sj associ6e £ une maille dSj. 



2847051 



14 



15 



Dans le cas de la propagation d'une onde electrique, les 
coefficients c if j, Vj et V(Mi), respectivement de la 
matrice d' interaction de la premiere et de la seconde 
matrice colonne, sont donn6s par : 



2n£o MjSj 



V(Mj) = Uj 

ou : 

- eo est une constante dielectrique, 

MjSj est une distance mesuree en valeur alg^brique, 

10 - qj correspond a une charge electrique caract^risant une 
source Sj, et 

- Ui correspond a un potentiel 61ectrique au point M±. 



Dans le cas de la propagation d'une onde magn§tique, 
1' expression de ces coefficients est la suivante : 



j =?j 
V<Mj) = 0j 



20 



oil : 

- ju 0 correspond & la permeabilite magnetique du milieu ou 
se situe le point Mi, 
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- q>j correspond au flux magnetique associe a la source 
0± correspond au potent iel magnetique au point Mi. 

Dans le cadre de la propagation d'une onde ultrasonore, 
ces coefficients sont donnas par : 



2tc M^j 



Vj=Vj 



[5] 



V(Mi) = Pi 

dans laquelle : 

- i 2 — 1. 

a> est la pulsation de l'onde sonore ; 

- p est la density du milieu dans lequel se situe le 
point Mi ; 

- le vecteur vj correspond a la vitesse sonore issue de 
la source Sj ; 

k correspond au vecteur d'onde de l'onde sonore ; et 

- Pi correspond & la pression acoustique generee par la 
propagation de I'onde ultrasonore au point Mj,'. 

Dans l 1 expression des coefficients ci,j, le terme dSj 
correspond k la surface de 1 ' echantillon associ§ a la 
source Sj, De preference, le maillage d'une surface au 
sens de 1'etape a) du proc^de selon 1' invention, est 
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choisi de sorte que chaque maille comporte une meme 
surface dS=dSi=dS2=~.=dSj . 

On remarque en particulier dans 1' expression des 
coefficients ci,j qu'ils dependent du produit scalaire 
entre le vecteur d'onde et le vecteur MjSj . Ainsi, pour 
des ondes ultrasonores, on tient compte d'un d§phasage 
entre les chemins qui lient chaque source Sj k un point de 
l'espace tridimensionnel M, ce d6phasage 6tant du a une 
difference de marche entre les rayons partant de chaque 
source et arrivant au point M (comme le montre la figure 
4A) . En particulier, on comprendra que 1' angle d' incidence 
d'un tel rayon est pris en compte dans 1' expression des, 
coefficients de la matrice d' interaction F. 

Bien entendu, dans le cadre de la propagation d'une onde 
electromagnetique de frequence elevee, done de longueur 
d'onde courte, qui differe du cadre 61ectrostatique ou 
magn6tostatique ci~avant, on tient compte du tefme de 
propagation exp(/£.r) dans 1' expression de la matrice 
d 1 interaction, par rapport & la geometrie du probl£me a 
resoudre, comme dans le cadre de la propagation d'une onde 
ultrasonore ci-avant (relation [5]). 

Ainsi, le syst^me matriciel de l 1 equation [1] permet 
d'estiraer, a partir d'une matrice d' interaction F et d'un 
vecteur comportant les valeurs Vj associ^es aux sources 
Sj, les coefficients d'un vecteur (matrice colonne) 
comportant les valeurs de grandeur physique V(Mi) aux 
points de I'espace Mi. 
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Pour determiner les valeurs des sources v jf on applique le 
systeme matriciel suivant : 



r V(P,) > 






V(P 2 ) 


= F x 




,V(P N ), 







[6] 



ou : 

- les coefficients de la matrice d' interaction F 
s'expriment par avec Cj-f{BSj) 

- les indices i et j correspondent respectivement £ la 
i** 6 ligne et la j 6me colonne de la matrice d' interaction 
F. Cette matrice d 1 interaction comporte, pour la 
determination des valeurs associees aux sources v-j, N 
lignes et N colonnes, en rappelant que N est le nombre 
total de mailles a la surface de 1' obstacle ; 

- les points Pi correspondent au sommet des hemispheres 
HEMi de la figure IB. 

La mise en oeuvre de l'6tape c) du procede au sens de la 
pr6sente invention correspond a calculer une condition aux 
limites pour les points Pi, de propri6t6s connues, comme 
on le verra plus loin. 

Le systeme matriciel de l f equation [6] devient alors : 
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hi 

v 2 


= F" 1 x 


V(P 2 ) 






,V(P N ), 



[7] 



ou : 

- F" 1 correspond a l 1 inverse de la matrice d' interaction 
5 F ; et 

- les valeurs V(P ± ) sont pr6d6termin§es, en fonction des 
conditions aux limites pr6citees. 

On determine ainsi les valeurs de source Vj. 

10 

A partir de l f estimation de ces valeurs de source Vj, on 
peut calculer la grandeur physique scalaire en un point M 
quelconque de l'espace tridimensionnel, & partir de la 
relation : 

15 

Pour obtenir cette expression de la grandeur scalaire 
V(M), la matrice d' interaction F peut ne comporter qu'une 
20 ligne de coefficients Cj, avec : 
Cj = fWSj) , 
mais comporte tou jours N colonnes. 

En.se r£f6rant a nouveau a la figure 1A, on comprendra que 
25 la surface de 1' obstacle OBS recevant l'onde qu'6met 
1' element rayonnant ER agit, elle-meme, comme une surface 
active re-§mettant une onde secondaire (par exemple par 
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reflexion) . Chaque source Si represente une contribution a 
1' emission de cette onde secondaire. 

En outre, pour tenir compte a la fois de la presence de 
I'onde principale et de la presence de I'onde secondaire 
au point M, on estime la contribution de I'onde principale 
et la contribution de I'onde secondaire au point M par le 
syst^me matriciel : 









V, > 


V(M) = Fx 




+ Fx 


v'2 











ou : 

- F» est la matrice d' interaction entre la surface de 
1 1 element rayonnant et le point M ; 

- v'j (avec j=l,2, 3,...,N' ) est la valeur des sources 
attributes a chaque 6chantillon de surface dS'j de 
l'616ment rayonnant, N 1 etant le nombre total de 
mailles choisi pour la surface active de 1' element 
rayonnant ER. 

Les coefficients de la matrice F' sont encore fonction de 
la distance MS'j, ou S'j sont les sources affectees a 
chaque tchantillon de surface dS'j de l'.616ment rayonnant. 

Selon une caract6ristique avantageuse, les valeurs des 
sources de 1" obstacle Vj sont determin6es en fonction des 
valeurs des sources de 1" element rayonnant v'j, lesquelles 
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sont elles-m^mes calcul6es comme on le verra plus loin en 
reference aux figures 4A, 4B, 5A et 5B. 

On se r6fere maintenant a la figure 2A, dans laquelle on 
5 affecte trois sources & chaque 6chantillon de surface dSi, 
en vue d'estimer une grandeur physique vectorielle V(M) , 
en un point M de l'espace tridimensionnel. 

En effet, on comprendra que pour estimer la grandeur 
10 vectorielle, par ses trois coordonn6es dans l'espace x, y 
et z, le nombre d' equations a rtsoudre par rapport au 
systdme matriciel precedent doit etre multiplie par trois. 
Ainsi, la matrice F" 1 de la relation [7] doit comporter 
trois fois plus de lignes que precederament . La matrice 
15 d f interaction F doit, elle-meme, comporter trois fois plus 
de colonnes que precedemment et, k cet effet, on pr6voit 
avantageuseraent trois sources par maille lorsqu'il s'agit 
de determiner les coordonn6es dans l'espace 

tridimensionnel d'un vecteur V(M) . 

20 

En se referant k la figure 2B, les trois sources SA if SB ± , 
SC ir attributes k un echantillon de surface dSi sont de 
positions respectives determinees comme on l'indique ci- 
apres. Telles que representees sur la figure 2B, les trois 
25 sources SAi, SB±, SCi sont coplanaires et le plan 
comportant ces trois sources comporte en outre la base de 
1' hemisphere HEMj.. L T hemisphere HEMi est construit comme 
indique ci-avant (de meme surface que la surface de la 
maille), avec toutefois le centre de l'hemisphere qui 
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correspond ici au barycentre des trois sources SA±, SB± et 
SC ± . 

On entend par "centre de 1 1 hemisphere" le centre du disque 
5 qui constitue la base de 1 'hemisphere • 

Les trois sources qui sont attributes ci 1 1 echantillon de 
surface dSi sont dispos6es aux sommets d'un triangle 
equilateral dont le barycentre G± correspond au centre de 
10 1 T hemisphere. De preference, chaque source SAi, SB ± et SC± 
est disposee sur le milieu d'un rayon R± de 1 'hemisphere, 
Ainsi, les droit es qui relient le barycentre G± a chaque 
source sont 6cartees angulairement de 120°. 

15 En se referant & la figure 2C, 1 ' orientation angulaire des 
triangles que forment les triplets de sources est choisie 
aieatoirement, d'un echantillon de surface a l'autre, 
Avantageusement, on 6vite ainsi des artefacts de sur- 
periodicite, dans 1' estimation de la grandeur vectorielle 

20 au point M, qui pourraient r6sulter du choix d'une meme 
orientation angulaire de ces triangles. 

En reference avec les differents types d'pndes indiques 
precedemment, la grandeur vectorielle V(M) & estimer peut 
25 etre : 

- un champ elect rique, dans le cadre de la propagation 
d'une onde 61ectrique ; 

un champ magnet ique, dans le cadre de la propagation 
d'une onde magnet ique ; et 
30 - une vitesse du son au point M, dans le cadre de la 
propagation d'ondes ultrasonores . 
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Pour determiner les valeurs associees a chaque source SAi, 
SBi, Sd, on met en forme le systeme matriciel selon la 
relation suivante : 



fv x (P,)l 

V X (P 2 ) 
V x (Pn) 

Vy(Pl) 
Vy(P 2 ) 

Vy(P N ) 

V z (Pi) 
V Z (P 2 ) 

,V Z (P N ) 



fvAi > 
VA2 

vA N 

vBj 

vB 2 

vB N 

vC 2 



[11] 



10 



On remarque, en particulier, que la matrice d' interaction 
Fy est de dimensions 3Nx3N, ou N est le nombre total 
d'echanti lions de surface. La matrice d' interaction 
s'exprime ici par la relation : 



Fv = 



N 



N 



N 



N)c£(iJ) C|(i,j) C£(i,j) 

N^(Uj) c^aj) cy(i,j) 

N^(i,j) C|(i,j) C^(i,j) 



[12] 
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Les coefficients de cette matrice s'expriment par : 

C v a iUj) = f u \d{Pi,Soj)] [131 

avec a = A, B, C 

i = 1, 2, N 
j - 1, 2, N 
u = x, y, z. 

et sont encore fonction d'une distance s^parant le point 
de contact Pi de l'une des sources SOj (<r = A, B ou C) 
d'un triplet associe a un 6chantillon de surface dS j . 

Par inversion de la matrice d' interaction Fy , on 
determine ainsi les valeurs vcrj associ<§es a chaque source 
S<Tj, en appliquant des conditions aux limites sur les 
valeurs du vecteur V aux points P±. Ces conditions aux 
limites imposent une valeur du vecteur V, selon ses trois 
coordonn£es V x (Pi), V y (Pi) et V 2 (Pi). 

Une fois que ces valeurs de sources vOj sont ainsi 
determinees, on calcule facilement l 1 expression du vecteur 
V en un point quelconque M de l'espace, par la relation : 

V(M) = V X (M) x + V y (M) y + V Z (M) z [14 ] 

V X {M) = ^f x [(d(M 9 S<7j)}.vcrj 
<T=A 9 B,C 
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&~A 9 B 9 C 

V Z {M)= %f z [(d(M 9 S*j)}.v*j 
a=A,B,C 

5 x, y et z correspondent a des vecteurs unitaires portes 
par les axes x, y et z de l'espace tridimensional. 

Ainsi, la matrice d' interaction F^ , lorsqu'elle est 
appliquee a un point quelconque M de I'espace, ne comporte 
10 finalement que trois lignes associ^es chacune & une 
coordonnSe de I'espace x, y ou z. 

Pour differents types d f ondes, les valeurs des sources vcr-j 
sont, comme prec^demment, une charge elect rique pour une 
15 onde 61ectrique, un flux magn6tique pour une onde 
magnetique, une vitesse du son pour une onde ultrasonore. 

Plus precisement, les coefficients de la matrice 
d' interaction Fy se d6terminent a partir des relations 
20 [3], [4] et [5] prec^dentes en pr<Scisant toutefois que : 

V(M) = -i^[V(M)] [15] 

V(M) 6tant la grandeur scalaire calcul6e pr6c6demment par 
25 l 1 equation [8] . 



25 



Ainsi, pour 1 ' estimation d'une grandeur vectorielle V au 
point M et pour les types d'onde cit6s ci-avant a titre 
d'exemple (electrique, magnet ique et ultrasonore) , les 
coefficients de la matrice d f interaction Fy sont 
inversement proportionnels au carr£ d'une distance 
s^parant chaque source du point M, alors que pour 
1 'estimation d'une grandeur scalaire V en un point M de 
l'espace, les coefficients de la matrice d' interaction F 
sont simplement inversement proportionnels & cette 
distance. Chaque distance implique l'une des sources d'un 
triplet d'un echantillon de surface et un point M de 
l'espace. La matrice d' interaction Fy comporte alors 3N 
colonnes lorsqu'il s'agit de prendre trois sources par 
echantillon de surface, alors que la matrice d 1 interaction 
F pour 1» estimation de la grandeur scalaire ne comportait 
que N colonnes puisque seule une source par echantillon de 
surface etait necessaire. 

Plus g6neralement, on attribue une source par echantillon 
lorsque 1' on connait des conditions aux limites pour une 
grandeur scalaire et trois sources par echantillon lorsque 
l'on connait plutdt des conditions aux limites pour une 
grandeur vectorielle. 

On se r^fere maintenant a la figure 3A pour d^crire, & 
titre illustratif , une application du proc6d6 selon 
1' invention a 1' estimation d'un potent iel 61ectrique en un 
point M de l'espace tridimensionnel, situ<§ entre deux 
armatures d'un condensateur . Les armatures de ce 
condensateur sont portees & des potentiels respectifs VI 
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et V2. La mise en ceuvre de l'etape a) consiste d'abord a 
mailler les surfaces respectives des deux armatures. Dans 
l'exemple represents* sur la figure 3A, on n'a repr6sente 
que deux mailles pour chaque armature, simplement a titre 
illustratif. 



10 



L' application de l'etape b) consiste a mettre en forme le 
systeme matriciel impliquant la matrice d' interaction F et 
le vecteur colonne comportant les valeurs des sources Si a 
S 4 . La multiplication de ces deux matrices permet 
d'obtenir un vecteur colonne comportant les valeurs du 
potentiel en un ou plusieurs points M de l'espace. 



La mise en ceuvre de l'etape c) du procede selon 
15 1' invention consiste a appliquer le syst&me matriciel aux 
points de contact des hemispheres Pi a P 4 , de chaque 



Schantillon de surface dSi a dS 4 . 
relation suivante : 



II en re suite la 



20 



v(p 2 ) 

V(P 3 ) 

Iv(p 4 )J 



2xe 0 



1 

P,Si 
1 

P 2 S, 

1 

P3S1 

1 

\ v P 4 S 1 



1 

P1S2 

1 

P2S2 
1 

P3S21 

1 

P4S2 



1 

P1S3 

1 

P2S3 

1 

P3S3 

1 

P4S3 P4S, 



P1S4 

1 

P2S4 

1 



P3S4 

1 



4*4 J 



v 2 
v 3 
^ v 4 



[16] 



avec V(Pi) = V(P 2 ) - Vi 

V(P 3 ) = V(P 4 ) = V 2 
Vi = qi, v 2 = <12, v 3 = q 3 , v 4 = q 4 
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Ici, la condition aux limites impose que la valeur du 
potentiel aux points de contact Pi et P 2 correspond au 
potentiel VI de la premiere armature. De m£me, le 
potentiel electrique aux points de contact P3 et P4 
5 correspond au potentiel Electrique de la seconde armature 
V2. Par inversion de la matrice d' interaction appliqu6e au 
point de contact Pi, on determine les valeurs des sources 
vi qui correspondent , comme exprime dans la relation [16], 
a des charges 61ectriques qi. 

10 

Les coefficients de la matrice d* interaction - — * ~ sont 

2^ 0 p i s j 

parfaitement connus, puisque les positions des sources Sj 
et les positions des points de contact Pi sont determinees 
au prealable, comme repr<§sent<§ sur la figure IB. 

15 

L 1 expression du potentiel electrique V(M) au point M entre 
les deux armatures est finalement donn6e par 
l 1 expression :. 

20 V (M) = ._Lf-SL + -Si-+-aa-+-a4»Y [i 7] 

Ittsq (MS] MS 2 MS 3 MS 4 J 

On se r6f 6re maintenant a la figure 3B sur laquelle on a 
repr6sent6 les m§mes armatures que sur la figure 3A, avec 
sensiblement un meme maillage, mais dans le but, ici, 
25 d'estimer une grandeur vectorielle correspondant au champ 

electrique E(M) , au point M de l'espace tridimensionnel . 
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On peut appliquer les relations [11] & [15] pour estimer 
la valeur du champ electrique au point M, avec, dans la 
relation [13] : 



1 



[18] 



avec 

tr = A, B, C 

u = x, y, z 

1 - 1, 2, 3, 4 

10 j = 1/ 2, 3, 4, 

Toutefois, les valeur s du champ 61ectrique au point de 
contact Pi restent a determiner dans la relation [11] . 

15 On introduit alors une loi g£n£rale predeterminee du 
comportement du champ (en reflexion, en transmission, ou 
autre) au niveau de la surface de 1' obstacle (des 
armatures dans 1' exemple du condensateur pr§cit6) , pour 
connaitre les valeurs des sources Va j . 

20 

Par exemple, si I'onde 61ectrique est totalement r6fl6chie 
par la surface d'un obstacle (par exemple l'une des deux 
armatures) , le champ electrique en un point de contact P± 
est normal a la surface dSi et ses composantes E x et E y 
25 sont nulles. A titre illustratif, si la surface de 
1' armature n'etait representee que par un seul 6chantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient toutes 6gales entre elles & une meroe 
valeur +q. 
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Au contraire, si le coefficient de reflexion est 
pratiquement nul h la surface dS ir la composante du champ 
electrique E 2 au point P± est nulle, ce qui correspond 

5 bien au cas ou le champ est sensiblement tangent £ la 
surface dSi. Ainsi, & titre illustratif, si la surface de 
1' armature n'etait representee que par un seul echantillon 
de surface a trois sources, les valeurs de ses sources vA, 
vB et vC seraient respectivement, par exemple, +q, +q et 

10 -2q. Par exemple, dans le cadre de la propagation d'une 
onde magnetique, si la surface d'un capteur a courants de 
Foucault (avec une composante normale du champ magnetique 
nulle) etait representee par un unique echantillon de 
surface, les flux magnetiques des trois sources associees 

15 a cet echantillon de surface seraient +9, +9 et -2<p. 

On comprend ainsi qu'avec les trois sources par 
echantillon dS±, on peut definir, par exemple en fonction 
de la ponderation de chaque source, une orientation 
20 quelconque du champ k la surface de 1' obstacle. 

Bien entendu, cette demarche suppose que le coefficient de 
reflexion R d'un obstacle soit connu auparavant . En 
particulier, il peut etre avantageux de comparer une 

25 simulation et une mesure exp6rimentale pour detecter, a la 
surface d'un obstacle, des inhomogeneites ou des impuretes 
qui correspondent a des points de la surface de cet 
obstacle qui ne v6rifient pas les valeurs du coefficient 
de reflexion R imposees & chaque point predetermine Pi de 

30 1' obstacle. 
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On peut affecter ainsi une valeur pr§determin6e du 
coefficient de reflexion a chaque point Pi de la surface 
de 1' obstacle. A cet effet, on introduit une matrice R qui 
est representative du coefficient de reflexion en chaque 

5 point P±. Pour une interaction entre un element rayonnant 
et un obstacle, on peut ainsi exprimer le systeme 
matriciel de la relation [9] autrement, c 1 est-£-dire en 
donnant une expression unique de toutes le? sources du 
systeme (a la fois de 1' obstacle et de 1' element 

10 rayonnant), comme indique ci-apr£s. 

Dans ce qui suit, on indique que : 
■ - F(P) est la matrice d ? interaction de 1' obstacle OBS 
appliquee aux points Pi de la surface de 1' obstacle 
15 OBS ; 

- F(P') est la matrice d' interaction de l f obstacle OBS 
appliqu6e aux points P'± de la surface de l v 616ment 
rayonnant ER; 

- F 1 (P) correspond & la matrice d' interaction de 
20 1* 616ment rayonnant ER appliquee aux points Pi de la 

surface de 1' obstacle OBS ; 

F' (P') correspond & la matrice d f interaction de 
l f Element rayonnant ER appliquee aux points P ± de la 
surface de l'61§ment rayonnant ER ; 
25 - v' correspond au vecteur colonne comportant les valeurs 
des sources S'i de 1' element rayonnant ER ; et 
v correspond au vecteur colonne comportant les. valeurs 
des sources Si de l f obstacle OBS. 

30 Sur 1' obstacle, la contribution de l'onde emise par 
l f element rayonnant ER s'exprime par : 
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V , (P) = F , (P).v f 



[19] 



La contribution de l'onde secondaire renvoyee par 
1* obstacle OBS s'exprime, par definition, par la 
relation : 



Or, dans l'exemple repr^sente sur la figure 4A, l'onde 
secondaire correspond simplement a une reflexion de l'onde 
principale. Ce qui s'exprime par la relation : 



ou R correspond & une matrice de reflexion dont chaque 
coefficient represente la contribution & 1' emission, par 
reflexion, de l'onde secondaire, par chaque source Si (ou 
Sai, dans le cadre d'une estimation d'une grandeur 
vectorielle) de 1' obstacle OBS, 

Des trois relations [19], [20] et [21] , on dSduit 
l 1 expression du vecteur colonne v comportant les valeurs 
des sources de 1' obstacle, & partir du vecteur colonne v 
comportant les valeurs des sources de l'616ment rayonnant, 
par la relation : 



V(P) = F(P).v 



[20] 



V(P) = R V'(P) 



[21] 



v = [F(P)]- ! .R .[F'(P)].v' 



[22] 
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Par ailleurs, pour une estimation fine des grandeurs 
scalaires ou vectorielles au point M, notairanent pour tenir 
compte de reflexions multiples, il est avantageux de tenir 
compte de la contribution du rayonnement par 1' obstacle, 

5 au niveau de la surface de 1' element rayonnant ER. A cet 
effet, on tient compte, dans l f estimation des conditions 
aux limites a la surface de l'element rayonnant ER (aux 
points P T i) de la contribution du rayonnement des sources 
S 1 i de 1 f element rayonnant et de la contribution de 

10 1' Emission de l'onde secondaire par les sources S± de 
1' obstacle, par la relation : 

VjCP'^FCP^v + FCPOv' [23] 

15 On peut ainsi ajuster, grace a la relation [23] , les 
valeurs de sources S'i de l f Element rayonnant ER f en 
tenant compte de la reflexion de 1' obstacle OBS, selon la 
relation suivante : 

20 V T (P f ) = \f(P') . \F(P) \ X . R . [F'(P) ] + F'(P') }.v' [10] 

Ainsi, on impose simplement des conditions aux limites 
pour 1' element rayonnant, pour en deduire les valeurs des 
sources v f i. En pratique, on procedera pr6f erentiellement 
25 comme suit : 

- apr6s maillage des surfaces, on determine la position 
des points Pi et P v i et des sources Si et S'i; 

- en fonction du type d'onde en jeu, on determine les 
coefficients des matrices P(P) , F' (P) , FfP 1 ) et F ' ( P 1 ) ; 
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- en fonction d'une loi de reflexion de l 1 obstacle, on 
determine les coefficients de la matrice de reflexion 
comme dans l'exemple donn6 plus loin pour une onde 
ultrasonore; 

5 - en fonction de conditions aux limites sur l'gl-gment 
rayonnant (dont le comportement est g§neralement connu 
pour un probleme donne) , on determine les valeurs du 
vecteur Vr{F) aux points P'i de la surface de 1' element 
rayonnant et on en deduit les valeurs des sources S'i de 

10 1' element rayonnant par inversion de la relation [10]; 

- on en deduit aussi les valeurs des sources Si de 
l v obstacle par application de la relation [22]; 

- une fois les valeurs de toutes les sources S L et S T i 
d6termin§es, on peut appliquer le syst^me matriciel donne 

15 par la relation [9] a tout point M de l'espace, en 
appliquant a ce point M les matrices d' interaction F et F 1 
(impliquant la position du point M et les positions des 
sources respectives S± et S'±) . 

20 En se r^ferant a nouveau a la figure 4A, on considere que 
1' obstacle OBS repr£sente simplement une interface entre 
deux milieux Ml et M2, formant ainsi un dioptre qui peut 
etre plan, tel que represent^ dans i T exemple de la figure 
4A, mais aussi incurv6 ou de forme generale quelconque. 

25 Les coefficients de reflexion R± associes £ chaque point 
Pi dependent, dans le cadre de la propagation d'une onde 
ultrasonore ou elect romagnetique de haute frequence, de 
l f angle d' incidence pi du rayon issu de la source Si, au 
point de l'espace tridimensionnel M. 



30 
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Pour une onde ultrasonore, 1' expression des coefficients 
de reflexion Ri est donn6e par : 



ou : 

- ci est la vitesse du son dans le milieu Mi ; 

- c 2 est la vitesse du son dans le milieu M 2 ; 

- px est la density du milieu Mi ; 
10 - p2 est la density du milieu M 2 . 

Dans cette expression [24] , le terme cosPi peut simplement 
etre estime en fonction des coordonnees dans I'espace du 
point M et du point repr6sentant la source Si- 

15 En se r§f£rant maintenant a la figure 4B, la meme 
estimation peut etre men6e pour un point M situ6 dans le 
milieu M 2 . Dans ce cas, l'onde que re<?oit le point M est 
une onde transmise par l'obstacle OBS- En particulier, on 
remarque que les sources de 1' element rayonnant ER ne sont 

20 plus actives, du fait de 1 1 occultation de 1' element 
rayonnant ER par l'obstacle OBS. En transmission, le 
raisonnement s' applique comme pr6c6demment avec une 
condition aux limites imposee aux points Pi par les 
valeurs des coefficients de transmission Ti associ6s a 

25 chaque point Pi- Dans le cadre de la propagation d'une 
onde ultrasonore, chaque coefficient de transmission Ti 
est donn£ par la relation : 
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2p2 c 2 COS A 



[25] 



1/2 



p 2 c 2 cos# +PiCj l-^T + ^i-cos 2 ^ 
L c l c l 



Comme indiqu6 pr6c6demment , les termes cosPi peuvent etre 



sources Si et du point M. 

Pour es timer les valeurs de sources S± de 1' obstacle OBS, 
on applique la relation [22] en remplagant toutefois la 
matrice de r6flexion R par la matrice de transmission T : 



Dans le cadre de la propagation d'une onde ultrasonore, 
les coefficients des matrices R et T s'estiment pour 
chaque source Si et pour chaque point Pi. En particulier, 
chaque coefficient T±,j ou R±,j de la matrice T ou de la 
matrice R (ou i correspond & la i feme ligne et j correspond 
a la j 6me colonne) s ! exprime en fonction d'un angle yf?ij 
entre une normale a la surface de 1' obstacle au point Pi 
et une droite passant par le point Pi et par une source 
Sj. On peut ecrire ainsi, de maniere gen6rale, les deux 
relations exprimant les valeurs des coefficients des 
matrices T et R par les relations respectives suivantes : 



determines en fonction des coordonn^es respectives des 



v = [F(P)] -1 T[F(P)]v' 



[26] 




ft (<x>s#j) 



[27] 
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R ifj = f r (cos/?ij) [28] 

ou ft est donne par la relation [25] et f r est donne par 
la relation [24] . 

De fagon plus gen6rale, en reference aux figures 4A et 4B, 
on indique que, si l'on considere 1' obstacle comme un 
materiau plein representant un milieu M2 distinct d'un 
milieu Ml dans lequel se propageait initialement l'onde 
principale : 

- pour une reflexion de l'onde principale sur l f obstacle 
en tant que milieu M2 (la surface de 1' obstacle formant 
dioptre entre les milieux Ml et M2) , les hemispheres HEM± 
sont orientes vers l f exterieur de 1' obstacle (figure 4A) ; 

pour une transmission de l'onde principale dans 
1' obstacle, les hemispheres HEMi sont orientes vers 
l'interieur de 1' obstacle (figure 4B) . 

On se refere maintenant a la figure 5A pour decrire le cas 
20 d'un 1' obstacle plan OBS de dimensions finies, excite par 
un element rayonnant ER, incline d'un angle predetermine 
par rapport a 1' obstacle OBS. Comme indique precedemment , 
pour une onde ultrasonore, on prendra en compte 
1 • inclinaison de 1' element rayonnant pour calculer la 
25 contribution de l'onde emise par 1' element rayonnant au 
point M. Par ailleurs, de fac?on particulierement 
avantageuse, on maille une surface qui englobe la surface 
de 1' obstacle (figure 5A) . Pour une tranche de l'espace 
delimit ee par 1 ' element rayonnant, d'une part, et 
30 l'obstacle, d'autre part (figure 5A) , on peut considerer 
trois types de sources : 
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- les sources S'i de 1' element rayonnant ER, 

- des sources SO if qui renvoient l'onde secondaire, par 
reflexion de 1' obstacle OBS, en fonction d f un certain 
coefficient de reflexion R de l 1 obstacle ; et 

• des sources SSi, qui ne renvoient pas d'onde secondaire 
et auxquelles un coefficient de reflexion nul peut etre 
attribue si 1' obstacle separe deux milieux d f indices 
identiques. Dans ce cas, ces sources SSi sont 
consid^rees comme "eteintes" dans la tranche de 
l'espace precit6e et ne sont pas prises en compte dans 
les calculs de la grandeur physique au point M de la 
figure 5A. En revanche, ces sources SSi peuvent etre 
actives par reflexion de l'onde principale si 
l'obstacle OBS separe deux milieux d' indices 
differents. 

Par ailleurs, pour estimer les grandeurs scalaires ou 
vectorielles associ^es a un point M d'un demi-espace 
delimite par la surface englobant l'obstacle OBS (a droite 
de la figure 5B) , on considere : 

les sources SO'i de l'obstacle, actives par transmission 

de l'onde principale, et 

les sources SSi, auxquelles on affecte maintenant un 
coefficient de transmission egal £ 1 si l'obstacle 
separe deux milieux de memes indices. Ces sources SSi se 
comportent finalement (aux angles d' incidence pres) 
comme les sources S'i de 1' Element rayonnant ER. 
Les sources S'i de 1' Element rayonnant peuvent alors §tre 
"eteintes" pour le calcul des grandeurs physiques dans ce 
demi-espace. 
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Pour calculer les valeurs v r des sources S'i de 1' element 
rayonnant (desquelles sont deduites les valeurs v des 
sources de 1' obstacle selon les relations [22] et [26] ), 
on appliquera simplement des conditions aux limites aux 
5 points de la surface active de 1' element rayonnant ER. Par 
exemple, pour une propagation d'onde ultrasonore, on peut 
indiquer que les vitesses sonores aux points de la surface 
de l'616ment rayonnant ER sont perpendiculaires a cette 
surface et de modules v 0 egaux entre eux. 

10 

De maniere g<§n6rale, on indique que l'espace 
tridirnensionnel peut etre ainsi d^coupe par des interfaces 
d^limitant des milieux de proprietes distinctes, chaque 
interface repr6sentant un obstacle au sens de la pr6sente 

15 invention. On calcule alors les grandeurs physiques dans 
chaque tranche de l'espace. Par exemple, dans le cadre de 
l'£tude d'une h£t6rostructure (£ plusieurs interfaces) on 
applique le proc6de ci-avant pour des tranches successives 
de l'espace en consid£rant deux interfaces : l'une 

20 representant un "Element rayonnant" au sens des figures 4A 
et 5A, par exemple par transmission d'une onde regue, et 
1' autre representant un obstacle recevant I'onde 
transmise. Avantageusement , on tient compte, pour chaque 
tranche de l'espace, des contributions de toutes les 

25 interfaces comme exprime par les relations [10] et [22].. 

La simulation de la figure 5C correspond, pour une onde 
ultrasonore, £ la situation des figures 5A et 5B en tenant 
compte : 

30 - de la contribution de 1' Emission de I'onde pri-ncipale 
par 1' Element rayonnant ER ; 
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- cie la contribution de la reflexion de cette onde 
principale par 1' obstacle ; et 

- de la contribution de la transmission de l'onde 
principale par 1' obstacle. 

5 

Les lignes de niveau de la figure 5C correspondent £ 
differents paliers de pression acoustique. L'el^ment 
rayonnant ER est dispose & 10 mm de 1" obstacle OBS et 
incline de 20° par rapport a ce dernier. On remarque en 

10 particulier des f ranges d f interferences dans une zone 
entre 1' obstacle OBS et 1 T element rayonnant ER. Une telle 
simulation peut avantageusement indiquer une position 
ideale d T un capteur ultrasonore. Ces capteurs ultrasonores 
comportent habituellement un tran'sducteur comme 616ment 

15 rayonnant actif et un detecteur pour mesurer les ondes 
ultrasonores regues. La simulation de la figure 5C peut 
ainsi indiquer en outre la forme ideale d'un capteur 
ultrasonore, selon les applications souhaitees, pour une 
forme d 1 obstacle donn6e. 

20 

La simulation de la figure 5C a 6t6 effectu6e grace a un 
calcul matriciel programme k l'aide du logiciel.de calcul 
MATLAB©. Le nombre de mailles total choisi pour 
1 'obstacle et pour l f Element rayonnant (ici, quelques 
25 centaines en tout) est alors optimist : 

- d'une part, pour limiter la dur<§e des calculs ; et 

- d 1 autre part, pour que la taille des mailles reste 
infer ieure a une demie longueur d'onde, de maniere a 
verifier le crit^re de Rayleigh. 
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On indique toutefois que, comme les elements k mailler 
dans la mise en ceuvre du proc6d6 selon l f invention sont 
simplement des surfaces, les temps de calcul sont bien 
moins longs que ceux necessaires dans la mise en oeuvre 
5 d'un precede de calcul de type "par 61&ments finis". 

La pr6sente invention peut ainsi se manif ester par la mise 
en oeuvre d'une succession d' instructions d'un produit 
programme d'ordinateur stocke dans la m^moire d'un disque 
10 dur ou d'un support amovible et se d6roulant comme suit : 

choix d'un pas de maillage notamment en fonction de la 

longueur d'onde de l'onde principale ; 

- determination des coordonn6es des sources S± et/ou S'i 
et des points de contact Pi et/ou P'i ; 

15 - choix d'un type d f onde en jeu et calcul des 
coefficients des matrices d' interaction appliquees aux 
points Pi et/ou P'i par la mise en oeuvre d'un logiciel 
de calcul matriciel ; 

- choix d'une loi de reflexion et/ou de transmission de 
20 la surface de l f obstacle et calcul des coefficients des 

matrices de reflexion et/ou de transmission ; 

- calcul des valeurs des sources S± et/ou S'i ; et 
calcul des grandeurs scalaires ou vectorielles en tout 
point de l'espace tridimensionnel. 

25 

A ce titre, la prfesente invention vise aussi un tel 
produit programme d 1 ordinateur . 

Bien entendu, la pr£sente invention ne se limite pas £ la 
30 forme de realisation d§crite ci-avant a titre d'exemple ; 
elle s'6tend a d'autres variantes. 
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Ainsi, on comprendra que, meme si dans les figures 
commences ci-avant on repr€sente k la fois la surface 
d'un obstacle et la surface d'un element rayonnant, la 

5 pr6sente invention s' applique aussi a l f estimation de 
grandeurs physiques dans le cadre d'une onde interagissant 
avec un obstacle et emise en champ lointain. Dans ce 
contexte, il n'est pas n^cessaire de mat^rialiser la 
surface d'un element rayonnant a mailler et les relations 

10 [8] et [14] ci-avant suffisent a determiner 1 ' interaction 
entre cette onde et 1' obstacle. 

On a indique ci-avant des Equations permettant de calculer 
les grandeurs scalaires ou vectorielles en un point M de 
15 l'espace, pour des ondes Electromagnet iques, ou encore 
acoust iques. Bien entendu, ces grandeurs peuvent etre 
estim^es pour d'autres types d f ondes , notamment pour des 
ondes thermiques, des ondes 61ectromagn£tiques impliquant 
des antennes radiof requences, ou autres. 

20 

Bien entendu, la pr6sente invention ne se limite pas a une 
application au contrdle non destructif , mais & tout type 
d f application, notamment en imagerie medicale, par exemple 
pour I 1 etude de microsystemes mettant en oeuvre une 
25 microscopie acoustique a miroirs mobiles . 

On a d6crit ci-avant des interactions entre une onde et un 
obstacle unique. Bien entendu, la presente invention 
s' applique a une interaction avec plusieurs obstacles. A 
30 cet effet, il convient simplement de mailler les surfaces 
de ces obstacles et d' additionner leur contribution pour 



2847051 



42 

1 1 estimation d'une grandeur vectorielle ou scalaire en un 
point quelconque de I'espace. De meme, comme indique ci- 
avant, la surface de 1' obstacle OBS peut etre plane, ou 
encore incurvee, ou encore de forme complexe quelconque. 

5 

Ainsi, dans le cadre d'une onde interagissant avec 
plusieurs obstacles dans l ! espace, une simulation 
equivalent e a celle representee sur la figure 5C 
permettrait de positionner des capteurs et/ou des Elements 
10 rayonnants en fonction de la configuration de ces 
obstacles , notamment pour une application a la 
determination de la position de hauts-parleurs dans un 
habitacle cloisonne, tel qu'un habitacle de vehicule 
automobile. 

15 

L'espace tridimensionnel peut etre d£coupe en une 
pluralite de regions, comme decrit ci-avant en reference 
aux figures 4A, 4B f 5A et 5B. Toutefois, pour qu'une 
surface de I'une desdites regions soit consideree comme un 
20 obstacle d'une onde principale, devenant actif par 
Emission d'une onde secondaire, 1 1 incidence de I'onde 
principale sur cette surface doit pref 6rentiellement 
rester inferieure ou 6gale £ 90°. 
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Revendi cations 

1. Proc6d<§ pour evaluer une grandeur physique associee a 
une interaction entre une onde et un obstacle , dans une 
region de l'espace tridimensionnel, dans lequel : 

a) on determine, par maillage, une plurality 
d' echantillons de surface (dS±) dont une partie au moins 
represente la surface d'un obstacle recevant une onde 
principale et emettant, en reponse, une onde secondaire, 
et l ! on attribue a chaque echantillon de surface au moins 
une source (Si) emettant une onde elementaire representant 
une contribution & ladite onde secondaire, 

b) on forme un syst&me matriciel comportant : 

- une matrice d' interaction (F(M)), inversible, 
appliquee a une region donn§e (M) de l'espace et 
comportant un nombre de colonnes correspondant & un 
nombre total de sources, 

une premiere matrice colonne dont chaque coefficient 
(vi) est associe a une source (S±) et caracterise 
1 ? onde elementaire qu'elle emet, 

- et une seconde matrice colonne, obtenue par une 
multiplication de la premiere matrice colonne par la 
matrice d' interaction et dont les coefficients sont 
des valeurs d'une grandeur physique (V(M)) 
representative de l'onde emise par 1' ensemble des 
sources en ladite region donnee <M) , 

c) pour estimer les coefficients de la premiere matrice 
colonne (v ± ) , on af fecte des valeurs de grandeur physique 
choisies (V(Pi)) a des points predetermines (Pi), associ6s 
chacun a un echantillon de surface (dSi) , lesdites valeurs 
choisies (V(Pi)) etant ranqees dans la seconde matrice 
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colonne, et l f on multiplie cette seconde matrice colonne 
par 1* inverse de la matrice d' interaction appliquee 
auxdits points predetermines (Pi), 

d) pour ^valuer ladite grandeur physique (V(M) ) 
5 representant l'onde emise par l f ensemble des sources dans 
une region donnee (M) de l'espace tridimensionnel, on 
applique la matrice d' interaction a ladite region donnee 
(M) et on multiplie cette matrice d' interaction par la 
premiere matrice colonne comportant les coefficients 
10 estimes a l'etape c) . 

2. Procede selon la revendication 1, dans lequel, pour 
evaluer une grandeur physique representative d'une 
interaction entre un element rayonnant une onde principale 
15 et un obstacle recevant cette onde principale , 

- A l f 6tape a), on determine en outre, par maillage, une 
pluralite d 1 fechantillons de surface (dS'i) representant 
ensemble une surface active de I 1 element rayonnant l'onde 
principale, et l'on attribue a chaque 6chantillon de la 

20 surface active au moins une source (S'i) emettant une onde 
elementaire representant une contribution a ladite onde 
principale, 

- on applique en outre les etapes b) , c) et d) aux 
echantillons de la surface active, et 

25 - on evalue ladite grandeur physique (V(M)) representant 
1 1 interaction entre l 1 element rayonnant et l 1 obstacle dans 
une r6gion donnee (M) de l'espace tridimensionnel, en 
prenant en compte la contribution, en ladite region donnee 
(M) , de l'onde principale emise par l 1 ensemble des sources 

30 de la surface active et la contribution de l'onde 
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secondaire emise par l 1 ensemble des sources de la surface, 
de 1 1 obstacle. 

3. Proc6de selon l'une des revendications 1 et 2, dans 
5 lequel chaque coefficient . de la matrice d f interaction, 
appliquee & une region donnee de l'espace, est 
repr^sentatif d'une interaction entre une source et ladite 
region donnee et la valeur de chaque coefficient est 
fonction d'une distance entre une source et ladite region 
10 donn6e. 



4. Proc6d6 selon l'une des revendications 1 a 3, dans 
lequel la matrice d' interaction appliquee, a l'etape c) , 
auxdits points predetermines (Pi) , comporte un nombre de 

15 lignes correspondant a un nombre total de points 
predetermines (Pi) . 

5. Proc6de selon l'une des revendications 1 a 4, dans 
lequel la grandeur physique k evaluer est une grandeur 

20 scalaire (V(Pi)) et, a l'etape a), on attribue une source 
unique a chaque 6chantillon de surface. 

6. Procede selon la revendication 5, dans lequel la 
matrice d 1 interaction ( F (M) ) appliquee, £ 1 1 6tape d) , & 

25 une region de l'espace (M) , comporte une ligne. 

7. Proced6 selon l'une des revendications 5 et 6, dans 
lequel chaque point predetermine (Pi) associ£ a un 
.6chantillon de surface (dSi) correspond a un point de 

30 contact entre cet echantillon de surface (dSi) et un 
hemisphere dont la surface est egale & la surface de cet 
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echantillon de surface, et de centre correspondant & une 
position de la source (Si) qui est attribute & cet 
Echantillon de surface. 

5 8. ProcEde selon 1'une des revendications 5 & 7, dans 
lequel : 

- l ! onde principale est une onde electrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge Electrique associEes chacune k une 

10 source, et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel electrique. 

9. Procede selon l'une des revendications 5 a 7, dans 
15 lequel : 

- l'onde principale est une onde magn§tique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnEtique associEes chacune a une source, 
et 

20 - les coefficients de la seconde matrice colonne sont des 
valeurs de potentiel magnetique. 

10. ProcEde selon l'une des revendications 5 k 7, dans 
lequel: 

25 - l'onde principale est une onde sonore, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune & une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
30 valeurs de pression acoustique. 
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11. ProcEdE selon l'une des revendications la 4, dans 
lequel la grandeur physique & Evaluer est une grandeur 
vectorielle (V(Pi)) exprimee par ses trois coordonnEes 
dans l'espace tridimensionnel, et l f on attribue, a 1 1 etape 

5 a), trois sources (SAi, SBi, SCi) a chaque Echantillon de 
surface (dSi) . 

12. Procede selon la revendication 11, dans lequel la 
matrice d' interaction <F V (M)) appliquEe, a l'etape d) , k 

10 une region de l'espace (M) , comporte une ligne pour chaque 
coordonnee (X,Y,Z) de I'espace. 

13. ProcEdE selon l'une des revendications 11 et 12, dans 
lequel : 

15 - les trois sources attributes a chaque Echantillon de 
surface sont sensiblement dans un meme plan et 
- chaque point predetermine (Pi) associE a un Echantillon 
de surface (dSi) correspond a un point de contact entre 
cet echantillon et un hemisphere dont la surface est Egale 

20 h la surface de cet Echantillon, et de centre 
correspondant a la position d'un barycentre des trois 
sources. 

14. ProcEdE selon la revendication 13, dans lequel les 
25 trois sources d T un meme echantillon de surface forment 

sensiblement un triangle Equilateral, et les triangles des 
echantillons de surface sont orientes sensiblement 
alEatoirement les uns par rapport aux autres. 

30 15. ProcEdE selon l'une des revendications 11 a 14, dans 
lequel : 
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- l'onde principale est une onde 61ectrique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de charge electrique associees chacune a une 
source, et 

5 - les coefficients de la seconde matrice colonne sont des 
valeurs de coordonnees d'un champ 61ectrique, 

16, Procede selon l'une des revendications 11 a 14, dans 
lequel : 

10 - l'onde principale est une onde magnetique, 

- les coefficients de la premiere matrice colonne sont des 
valeurs de flux magnetique associees chacune & une source,, 
et 

- les coefficients de la seconde matrice colonne sont des 
15 valeurs de coordonnees d'un champ magnetique. 

17, Proced6 selon l'une des revendications 11 a 14, dans 
lequel : 

- l'onde principale est une onde sonore, 

20 - les coefficients de la premiere matrice colonne sont des 
valeurs de vitesse de son associees chacune a une source, 
et 

- les coefficients de la seconde matrice colonne sont des 
valeurs de coordonnees d'une vitesse acoustique. 

25 

18, Proc<§d6 selon l'une des revendications pr6c^dentes, 
dans lequel, pour estimer la contribution de l'onde 
secondaire en ladite region donnee a l'etape d) , lesdites 
valeurs de grandeur physique (V(Pi)) choisies a l'etape c) 

30 sont fonction d'un coefficient predetermine de reflexion 
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et/ou de transmission de l ! onde principale par chaque 
echantillon de surface de I 1 obstacle. 

19. Procede selon la revendication 18 , prise en 
combinaison avec I'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond a une reflexion de 
l'onde principale sur 1' obstacle et l 1 hemisphere est 
oriente vers 1 ' ext£rieur de 1' obstacle. 

20. Proc6d6 selon la revendication 18, prise en 
combinaison avec l'une des revendications 6 et 12, dans 
lequel l'onde secondaire correspond h une transmission de 
l'onde principale dans I 1 obstacle et 1' hemisphere est 
oriente vers 1' inter ieur de 1' obstacle. 

21. Procede selon l'une des revendications 19 et 20, dans 
lequel, a l'etape c) , on determine les valeurs (v'i) 
associees aux sources (S'i) de 1' element rayonnant (ER) et 
on met en forme au moins : 

- une premiere matrice d' interaction (F(P)) repr^sentant 
la contribution des sources de 1' obstacle aux points 
predetermines de la surface de I'obstacle (Pi), 

- une seconde matrice d' interaction (F' (P) ) repr£sentant 
la contribution des sources de l'£16ment rayonnant aux 
points predetermines de la surface de I'obstacle (Pi), 

- une matrice de reflexion (R) ou de transmission (T) , 
dont les coefficients repr£sentent des coefficients de 
reflexion ou de transmission en chaque point predetermine 
(Pi) de I'obstacle, 

pour determiner les valeurs des sources de I'obstacle (vi) 
en fonction des valeurs des sources de 1' element rayonnant 
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(v'i) et d'une multiplication des premiere et seconde 
matrices d 1 interaction et de la matrice de reflexion ou de 
transmission. 

5 22. Procede selon la revindication 21, dans lequel, £ 
l'etape c) , on determine les valeurs (v'i) associees aux 
sources (S'i) de 1' element rayonnant (ER) en tenant compte 
de la reception de I'onde secondaire par 1" element 
rayonnant (ER) et en mettant en forme en outre : 

10 - une troisieme matrice d' interaction (F(P')) representant 
la contribution des sources de 1' obstacle aux points 
predetermines de la surface de l'felement rayonnant (P'i), 

et une quatrieme matrice d' interaction (F' (P f ) ) 
representant la contribution des sources de 1' element 

15 rayonnant aux points predetermines de la surface de 
1' element rayonnant (P'i) . 

23. Procede selon l'une des revendications 19 a 22 , dans 
lequel la surface de 1' obstacle correspond a une interface 

20 entre deux milieux distincts d'une heterostructure . 

24. Procede selon l'une des revendications precedentes, 
dans lequel 1 1 onde principale est une onde sonore et les 
coefficients de la matrice d' interaction sont chacun 

25 fonction d'un angle d' incidence d'une onde eiementaire 
issue d'une source en ladite region donnee (M) . 

25. Procede selon l'une des revendications pr6cedentes, 
dans lequel 1 ' onde principale est une onde sonore et, a 

30 l'etape a) , on choisit un nombre total d ' echantillons de 
surface (dSi) sensiblement en fonction d' une longueur 
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d'onde de l'onde sonore pour verifier le crit^re de 
Rayleigh. 

26. Proc§d6 selon 1 1 une des revendications pr6c§dentes, 
dans lequel on compare une pluralite de valeurs de la 
grandeur physique estimee a l'etape d) , obtenues pour une 
pluralite de regions de l'espace, pour selectionner une 
region candidate pour . la disposition d'un element 
rayonnant destine a interagir avec l 1 obstacle. 

27. Procede selon l'une des revendications 2 k 26, dans 
lequel 1' element rayonnant est un capteur, pour un 
controle non destructif, destin§ a analyser un objet 
formant un obstacle de l'onde principale. 
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